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Abstract 


In this paper we study the flow of an inviscid fluid composed by three different phases. The 
model is a simple hyperbolic system of three conservation laws, in Lagrangian coordinates, 
where the phase interfaces are stationary. Our main result concerns the global existence of 
weak entropic solutions to the initial-value problem for large initial data. 
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1 Introduction 

The theory of hyperbolic systems of conservation laws in one spatial dimension has reached in 
the last years a rather satisfactory level of completeness, as the reference book of Dafermos [12] 
witnesses. Among the several important results that have been proved, probably the greatest 
achievement concerns the global existence in time of weak solutions to the initial-value problem, 
as well as their uniqueness, continuity with respect to the data and viscous approximations. 
However, such results hold, in general, only for small initial data: the case of large data has 
been given no general and satisfactory answer. This paper focuses precisely on this issue in 
the case of a simple but physically meaningful system of three equations, for which we provide 
explicit conditions on the initial data in order to have global solutions. 

The system under consideration arises in the modeling of phase transitions for an inviscid 
fluid and is deduced by [T3|. If we denote by n > 0 the specific volume of the fluid, u the velocity, 
p the pressure and A the mass-density fraction of the vapor, the system is written as 


/ 


Vt - Ux 

< Ut+p{v,X}x 

At 

V 


= 0 , 
= 0 , 
= 0 . 


( 1 . 1 ) 
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As usual, here t > 0 denotes the time and x G M. The phase states of the flnid are modeled by 
the variable A, which ranges from 0 (pure liquid) to 1 (pure vapor) and allows for intermediate 
values in the interval ]0,1[ representing mixtures of the two pure phases. The model incorporates 
the state variable A in the pressure, which is dehned by 


p(u. A) = 


a2(A) 


( 1 . 2 ) 


where a(A) > 0 and is a function on [0,1]. We denote hy U = {v,u, A) the state variables 


and by O =]0,+cx)[x]R x [0,1] the state space. System (1.1) is strictly hyperbolic in with 
eigenvalues ei = ^2 = 0 , 63 = the hrst and the third characteristic helds are 

genuinely nonlinear, while the second one is linearly degenerate. 


The first result on the existence of global solutions to system (1.1), provided with suitably 


large initial data, is given in [3]; a different proof is given in [7]. In particnlar, in the case where 
A is constant, the classical result by Nishida [16] is recovered. 

The analysis of |3| is pursued and refined in pQ for initial data with A of Riemann type: 


A(x, 0) is constant for x / 0 with a jnmp at 0. In this case system (l.I) decouples for any t > 0 
into two p-systems connected by a phase interface at x = 0, because the discontinuities of A do 
not propagate. The special form of A allowed us to analyze in detail the effect of the nonlinear 
interaction of pressure waves through the phase interface, leading to refined sufficient conditions 
on the initial data for which solntions exist globally in time. 

A survey on couplings of two systems of conservation laws, with a focus on numerical ap¬ 


proximations, is given in mi; however, we emphasize that the above coupling for system ( | 1 . 1 [ ) 
is the physical coupling, where the interface is a contact discontinuity. We refer to [H [1] for 
further references on related results. 


In this paper we continue the analysis of system (1.1) by considering the case where the 


initial datum for A is piecewise constant with two jumps. Let the initial data be of the form 


Uo{x) = (uo(x),?Xo(x), Ao(x)) , 


for Ao(x) 


/ 

Xi if X < fl, 

< Xm if a < X < ^ , 
Xr if X > 6 , 


(1.3) 


where x G M and A^, Xm, Xr G [0,1] are constant. We define ag = a{X£), am = a{Xm), cir = a{Xr) 
and focus on the case 

am <xain{ai,ar] . (1.4) 


To give a flavor of the physical meaning of the problem, assume that a (A) is increasing (which 
is the physically meaningful case) so that (1.4) implies Xm < minjA^, A^}. Here, we are dealing 
with a one-dimensional fluid consisting of three homogeneous mixtures of liquid and vapor; the 
mixture in the region ]a, ^[ is more liquid than in the surrounding ones. This inclndes the case of 
a liqnid drop in a gaseous environment. The other cases am > max{af, a,.} (a bubble surrounded 
by liquid) and ae < am < a-r (or ai > am > o,r) are considered in a forthcoming work [2j. 

A similar model is studied in [TO]. There, the basic system (in Eulerian coordinates) has only 
two equations but is augmented with kinetic conditions deduced by the mass and momentum 
conservation at the interfaces, which make that model essentially equivalent to (1.1). However, 
the resnlts of m concern a general pressure law but small initial data. We refer also to mi 


2 













where an analogous system (in Lagrangian coordinates) is studied, in which the pressure in the 
region [a, 6] is a linear function of v. 

Notice also that (1.1), (1.3) can be interpreted as a perturbation problem of the steady 
solution given by the two parallel contact discontinuities located at x = a and x = 5, respectively. 
We refer to [Eld for the analysis of the perturbation of a single contact discontinuity. 

The main result of this paper is Theorem 2.1, that provides a wide class of large initial data 
for which the solution to the initial-value problem (1.1), (1.3) exists globally in time. Roughly 
speaking, the conditions on the data require that the total variations of po = p{vo, Ao) and Uq 
do not exceed a certain threshold depending on the sizes of the interfaces, see (2.3): the larger 
are the interfaces, the smaller must be the variations and conversely. Also, if the variations are 
sufficiently small then any size of the interfaces is allowed, provided that the stability condition 
(2.1) (that was missing in pQ) holds. Such a result was proved, to the best of our knowledge, in 
no related paper. Moreover, we point out that our conditions on the initial data are sufficiently 
flexible to allow the control of the variations in either of the three phases. 

For the proof of Theorem 2.1, two novel ideas are employed. The hrst one is a simplification 
in the dehnition of the functional F used to control the total variation of the solutions (see (5.1) 
and Remark 5.1), in which some nonlinear terms are dropped, thanks to a more careful use of 
nonlinear interactions involving phase waves. 

The second one is an original variant of the front-tracking algorithm [9], that is needed in 
order to ensure that the functional F is decreasing. Indeed, the classical front-tracking scheme 
prescribes two ways of solving the Riemann problem arising at an interaction: either by means 
of an Accurate solver or by a Simplified solver that exploits non-physical waves, which are used 
to prevent the possible blow-up in hnite time of the number of fronts and interactions. Here 
we provide an original definition of the Simplified solver, suitably designed for this problem. At 
any interaction of a small wave with an interface, the Simplified solver introduces stationary, 
non-entropic waves (associated to the integral curves of 0), which are formally computed 
as reflected waves. These waves “travel” with zero speed and then remain attached to the 
phase wave, thus forming a “composite wave”. Such a Riemann solver somewhat reminds of the 
famous Osher solver frequently used in Numerical Analysis, see |18L §12]. The idea of introducing 
stationary composite waves for the Simplified solver is also exploited in [T] where, however, the 
jump across non-physical waves is defined as in [9]. 

Notice that when A is constant, system (1.1) reduces to a 2 x 2 system and one can avoid the 
use of the Simplified solver (see mm- We point out that non-physical waves are also avoided 
in mm, but for different reasons: in uni, due to a particular solver and to the smallness of 
the data, while in m because of the assumption of linear pressure in the region [a, 5]. 

The paper is organized as follows. The main result is stated in Section In Section we 
first introduce four pre-Riemann solvers: one of them is used to define the composite wave, 


the other three are exploited either in the Accurate or in the Simplified solver. Proposition 3.7 


gives a unified approach to both the Accurate and the Simplified solver. Approximate solutions 
are defined in Section [4[ In Section [5] we introduce the main functional F and show that it 
is decreasing in time. Section deals with the convergence and consistency of the algorithm, 
together with a decay property of the reflected waves; we provide also a comparison with jlj 
which shows how the preceding result is improved. Finally, in Appendix [^ it is proved an 
alternative estimate concerning certain interactions solved by the Simplified solver. 
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2 Main Result 


Throughout this paper we assume (1.4) and call tj, the strengths of the two 2-waves, as in [3j: 


T] = 2 


dm 


C = 2 


CLr' O'Ti 


dr + dr, 


By (1.4) and for a£, dm, CLr in M"*" =]0, +oo[, one easily finds that 

77 <0, C>0, |r7|, Id G [0,2[. 

In what follows, we need that rj and ^ satisfy the stdbility condition 


( 2 . 1 ) 


When one of the two waves t] ox C vanishes, for example C = 0, then reduces to \r]\ < 2, 
which is always satisfied. The inequality (2.1) identifies a set 2? C [0,2[x[0, 2[ (see Figure [^, 
where we define a non-negative and continuous function T-L by 


«(|.,|,|C|) max| j + |'■d2)|,|/2■ 1 - (1 + |’!j|/2)|C|/2} 


( 2 . 2 ) 





Figure 1: The domain D in the (I 77 I, |d)-plane. 


Notice that 22 = 0 only when 77 = (^ = 0; it holds 22 (|r 7 |, 0 ) = \r]\ and 22(0, |d) = |CI- 
Moreover, we have that 22 (|t 7|, |d) tends to -|-oo when (|r 7 |, |d) tends to the curved edges of V. 
Following the notation in Figure we set 


£ = {(x, t) : X < a} , M. = {(x, 2) : fl < X < ^} , TZ = {(x, t) \ x > 6} . 


We denote po{x) = p {vo{x), Ao(x)) and TY{f,g) = TV f + TY g, for any / = f{x),g = g{x). 


Theorem 2.1. Assume and consider initial data ( |1.3[ ) with Vo{x) > v > 0, 

constant V. Assume also (1.4) and (2.1). There exists a strictly decreasing function 
for r > 0, with 


for some 
K, defined 


lim /C(r) = -1-00 , lim /C(r) = 0 , 

r—). 0 -|- 1—>--1-00 
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Figure 2: The regions £, TZ in the (x, t)-plane. 


such that if it holds 


TV 

x<a 


log(Po), - 
ae 


1 


+ 1 I I Th flog(po), — ) + TV ("log(po), —) < V ('H(|?7|, ICD) , 

1 + n[\r]\,\Q) a<x<6 \ QraJ Or J 

(2.3) 

then the Cauchy problem (1.1), (1.3) hasaweak entropic solution {v,u,X) definedfort G [0,+oo[. 
If I] = C = ^ same conelusion holds with K, (T-i{\rj\, |^|)) replaced by +oo in (2.3). 

Moreover, the solution is valued in a compact set and {v{-,t),u{-,t)) G L°°([0, oo[; BV(M)). 


The definition of the threshold function K, is given in ( |6.18 ) and is the same as [H (6.20)]. 

, ICI can be taken, the smaller the 
r/l, ICI, the larger can be the total 


Hypothesis (2.3) can be interpreted as follows: the larger 


total variation of po, Uo must be; vice versa, the smaller are 
variation of Po,Uo- 

Consider, now, the case when one of the two phase-waves tends to zero, say 

n\vi\(\) ~ 


0. Then 


TV 

x< 


H(0, Id) = Id and ( |2.3[ ) becomes formally 
1 


/(log(pd,-) + 


TV 


log(Po), —) + TV (log{po), < /C(|d), (2.4) 

am) V Or) 


1 -|- Id a<x<6 

which improves lU (2.3)] by allowing to take larger total variation of the data for x G]fl, ^[. 
Indeed, when It^I =0, we will see in Remark 6.4 that hypothesis (2.4) can be improved by 

1 


TV 


(log(po), 


— ) + TV 


log(Po), — 

Qj'p 


<m\), 


(2.5) 


1 + Id V 

by which the total variation can be taken larger on the entire interval ] — oo, ^[. 

Theorem 2.1 improves also the main result in [3| when restricted to the case of two contact 
discontinuities, not only because K, is sharper than // of [5l Theorem 3.1] (see Sectionbelow), 
but also because the total variation of the initial data (thanks to the coefficient of the middle 
term in (2.3)) can be larger in Mi. Recall that Mi is the more liquid region if a(A) is increasing. 
The asymmetrical character of (2.3) is due to the particular choice of the decreasing functional 
F used to estimate the total variation of the approximate solutions, see Section 

We also notice that a slight improvement of condition (2.3) in Theorem |2.1| would follow 


from the use of the Riemann coordinates, see Remark 6.3 


We conclude this section by extracting some more information from (2.3); with this aim we 
introduce the sub-level sets of F., 


Vr = 


,\C\)eV: ?^(|d,|d) <c} , OO, 
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see Figure]^ Since K, is decreasing, then for every {\ri\, |C|) G Vc condition (2.3) holds if 


TV 

x<a 



+ TV 

a<x<B 



+ TV 

X>l} 



< JC{c). 


In particular, we have /C(2) = 21og(2 + \/3)/3 and the domain 7)2 includes the segments [0,2[ 
on each axis. Therefore, for r/ = 0 or = 0 we recover a slightly better condition than [U (2.5)]. 
We notice that the 2-level set of Ti has a particular simple expression; it is the graph of the 
function C(|??|) = 2(2 - |r/|)/(2 |ry|). 



Figure 3: Sets of level c of the function V.: cases c = 1, 2, 3. 


3 The Riemann Problem 


In this section we collect some basic facts about system (1.1); we refer to dElIl! for more 
details and to PUS] for generalities on Riemann problems. As anticipated in the Introduction, 
in addition to the usual Lax waves used in the theory of conservation laws, we introduce suitable 
composite waves which sum up the effects of each contact discontinuity and of certain reflected 
waves. Finally, we present two Riemann solvers that use such composite waves. 

For i = 1,3, the f-th right shock-rarefaction curves <I>i through the point Uo = {vo, Uo, Ao) G 12 


for (1.1) are as in [1] 

V 1 -^ ^i{si){Uo) = {v, Uo + 2a{Xo)h{si), Ao) , v > 0 , i = 1,3 , 
where the strength gj of an f-wave is defined as 








and the function h is defined by 


(3.1) 


(3.2) 


h{e) = 


if e > 0 , 
sinh e if e < 0 . 


(3.3) 
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Rarefaction waves have positive strengths and shock waves have negative strengths. The z-th 
integral curve through ?7o G n is denoted by Ii{e){Uo), for e G M and z = 1, 3; two states U and 
Ii{e){U) are connected by an z-rarefaction wave iff e > 0. The wave curve corresponding to the 
second characteristic field through C/o G O is given by 






A G [0,1], 


and the strength of a 2-wave is 

_ g «(A) -a(Ao) 

' a(A)+a(Ao)- 


For starters, we prove a result similar to [H Proposition 3.2], For \± G [0,1], we use the 
notation a± = a(A±), p± = p{v±, A±). 

Proposition 3.1. Fix two functions 01,^3 that can he either the identity Id or the function h 
defined in (3.3). For any pair of states U- = (z;_, zz_, A_), C/+ = (z;+,n+,A+) G It, there exist 
unique ei, £3 G M such that: 


£3 - ei = ^ log 


( P± 

\P- 


o_ 6 'i(ei) + a+ 6 » 3 (e 3 ) = 


ZZ+ — U- 


(3.4) 


Proof. Let us call log(p+/p_)/2 =: A and {u+ — zz_)/2 =: B, since they are two constant 
quantities once we fixed U- and C/+. Thus, we have four possible cases to examine: 


£3 - ei = A , 
a-h{ei) + a+h{ez) = B , 

£3-Si= A, 
a-h{ei) + 0+^3 = B , 


£3-£i = A, 

a-£i + a +£3 = B , 

£3-£1= A, 

a-£i + a+h{£ 3 ) = B . 


(3.5) 

(3.6) 


System ( |3.5[ )i {9i = 63 = h) has already been solved in jH Proposition 3.2], while system ( 3 . 5)2 
is linear. As for (3.6), it suffices to study just one of the two systems, for example (3. 6)1 (the 
other one is analogous). In this case, setting k = a+/o_, it holds h{£i) + k£i = B/{a-) — kA. 
Thus, if G{x) := kx + h{x), we have G(ei) = B/{a-) — kA. Since G is invertible and onto, this 
gives £i = G~^{B/{a-) — kA). □ 


Remark 3.2. Notice that only system (3.5 )i always gives an actual Lax solution to the Riemann 
problem for (1.1) with initial data 

U{x,0) = 


U- if x < 0 , 
[/+ if x > 0 , 


(3.7) 


as the juxtaposition of a 1-wave of strength £i, a 2-wave 5 = 2(a+— a_)/(a+-|-a_) and a 3-wave 
of strength £ 3 , see Q Proposition 3.2]. In general, this is not true for the other three cases. 

When solving an interaction with a 2-wave 6 , we sometimes make use of a Riemann solver 
that attaches certain reflected waves to 5; the outcome is a stationary composite wave, which is 
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made of the composition of a wave related to an integral curve for the first characteristic field, 
the 2-wave 6 , and, finally, a wave related to an integral curve for the third characteristic field. 

We use the symbols ‘L’ to refer to the Lax curves and ‘I’ to refer to the Integral curves 
li, i = 1,3. Then, Proposition |3.1| allows us to give the following important definition. 


Definition 3.3 (Pre-Riemann solver). For any choice 0 / 61,63 as in Proposition 3.1. the Pre- 
Riemann solver : If x Q— t-Mx] — 2,2[x]R is the map defined by 


ReidsiU-, U+) — £ 3 ), 


(3.8) 


where £i ,£3 are as in (3.4) and 6 = 2(a+ — a_)/(a+ + a_). The two subscripts in 61,63 stand 
for the choice of 1,3-wave curves (L or I) along which £1,63 are taken. More precisely, it holds 
6 i = h for Si along Lax curves, while 6 i = Id for Si along integral curves. Then, by Proposition \3.1\ 
we get four Pre-Riemann solvers that we denote by Ru,, Rn, Rn o,nd Rjl, respectively. 


Notice that Ru, is an actual Riemann solver by Remark |3.2| and Rn, Rn are used in 

We do not assign any 


connection with the Simplihed Riemann Solver, see Proposition 3.7 


speed to £i ,£3 when they are taken along integral curves /; indeed, these waves shall be sticked 
to the phase wave and can be thought as being stationary. In particular, Ru is used to dehne 
composite waves as in the following dehnition. 


Definition 3.4 (Composite wave). A composite wave So = (dQ,S,dQ) associated to a 2-wave 6 
and connecting two states U- = (u_,u_, A_) and 11+ = (u+,u+,A+) of Ft, with A_ / A+, is the 
wave with zero speed defined by So = Ru{U-, U+). We write |(5o| = |(^ol + I'^ol- 

Notice that 5o reduces to a 2-wave as long as (5 q = (5g = 0. We denote by rjo and Co the two 
composite waves associated to rj and C, respectively; see Figure]^ 



Figure 4; The composite waves in the (x, t) plane: in (a) rjo and 
Co are drawn as three parallel close lines, while (b) is the auxiliary 
picture that is used to determine the states in the interactions, see 
Figure 


Remark that in Figure]^ (6) the r?o,Co components, i = 1,3, may be non-entropic waves: they 
are depicted as fronts with negative speed (i = 1) and positive speed (i = 3) in order to easily 
understand how to handle the interactions. 

In this way, we are left to deal with waves of family 1,3 and two distinct composite waves 
belonging to a fictitious 0-family. Notice also that, once we fix U-, the set of states U+ that can 














be connected to U- by a composite wave does not describe a curve in the {v, u) plane, but the 
whole half-plane u > 0. 

Before proceeding with the detailed description of the two new Riemann solvers, we insert 
here the following useful lemma. Now and then we will make an inappropriate use of the term 
‘waves’ to indicate both actual physical waves (i.e. connecting states that lie on a Lax curve) 
and not (i.e. when referring to states that lie on a general integral curve or on a combination of 
Lax curves and integral curves). 

Lemma 3.5 (Commutation of Lwaves). Let i = 1, 3 and ai, (3i G M. If two states U-, C/+ G H 
in the same phase (X- = A+J are connected by an i-wave of strength ai followed by an i-wave 
of strength Pi, then they can be connected also by an i-wave Pi followed by an i-wave ai. 

Proof. Assume i = 3 (the other case is analogous) and fix 6 ^^ , 9^^ to be either h or Id, see 
Figure]^ If U* = {v* ,u*,X*) is the final state reached starting from U- and moving first along 
P^ and then along 03 , then trivially it holds A* = A+ and 

V* = V- exp(- 2/33 - 203 ) = , u* =U--[- 2a_ , 

that means U* = U+. □ 



Figure 5: The commutation of i-waves: case i = 3, a 3 ,P 3 < 0, 
9^^ = h and 9^^ = Id. Here Um and Uq are the states connected 
to U- along the 3-Lax curve by 03 and, respectively, along the 3- 
integral curve by P^. 


Remark 3.6. When 


= h, Lemma 


3.5 


is a consequence of the invariance by translation 
of Lax curves for the p-system with 7 = 1; see m- 


Now, we are ready to describe the two Riemann solvers that will be needed in case of 
interactions with 770 and Co at positive times: we use an Accurate solver when the interacting 
wave has size bigger than a threshold p to be determined and a Simplified solver otherwise. 

We denote by 6 i (and ep the interacting waves (the waves produced by the interaction, 
respectively), for i = 1,3; note that, taking for simplicity 5 to be equal either to p or to 
we use the same notation 60 = ((5q,(5, <5q) (and Eq = (£o)<^)£o)) to denote both pQ and Co as 
interacting waves (and as outgoing waves, respectively). 
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Proposition 3.7. Let z = 1, 3 and consider the interaction at a time t > 0 of a composite wave 
So = ((5g,(5, 5 q) with an i-wave of strength Si; we refer to Then, the emerging Riemann 

problem with initial states U-,U+ can be solved by means of in one of the two following 
ways. Denote U- = Ii{Sq){U-) and U+ = Io{—Sq){U+). 

1. Accurate Riemann solver. The solution is formed by waves ei,eo,£'i, where (ei,<5,£ 3 ) = 
Rll{U-,U+) and eo = Sq. 


2. Simplified Riemann solver. We distinguish case i = 1 and i = 3.' 

i) for i = 1, the solution is formed by waves £i,£o such that (£i,<5,£ 3 ) = R]j{U-,U+) 
and £0 = (<^ 0 ,( 5 , 5o +£ 3 ); 

a) for i = 3, the solution is formed by waves £o ,£3 such that (£i,(i, £ 3 ) = Rjl{U-,U+) 
and £0 = (^0 + ^1) < 50 ) • 


In general, it holds 9i = h in all cases; for any Oj, j = 1, 3, j / i, chosen between Id and h, 
the following relations are verified: 


£3 - £1 


-5i ifi = l, 

S 3 if i = 3, 


a_6»i(£i) + a+6'3(£3) 


Moreover, in all cases the signs of 81,83 satisfy: 


a+9i{Si) 

a- 03 ( 53 ) 


if i = l, 
if i = 3. 


(3.9) 


sgn£j = sgn5i, sgn£j 


sgn S -sgn Si if i = 1, 
—sgn S ■ sgn Si if i = 3 . 


(3.10) 


Proof. In the interaction of an z-wave Si with a composite wave Sq, we look at the interaction of 
Si with the S component of Jq: indeed, Si crosses Sq, j = 1,3, j / z, without changing strength 
by m Lemma 5.4]. Then, we solve the Riemann problem with initial states by means 

of ^ 6 * 103 ) with 01,03 either Id or h. We proceed as follows. 


1. Accurate Riemann solver. After computing Rll{U-,U+) = (£i,(5,£ 3 ), we let £1 and £3 
commute with (ig and (5g respectively, in the sense of Lemma |3.5[ In this way, they are 
free to propagate as outgoing waves of family 1,3; see Figure]^ (a), ( 6 ) for a picture of case 
i = 3. Then, the resulting composite wave connects Up to Uq, where Up = ‘hi(£i)(17_) and 
Uq = $3(-£3)(t/+). Hence, £0 = Rii{Up,Uq) = {SiS,Sl) = Sq. 


2. Simplified Riemann solver. We have to distinguish between case z = 1 and z = 3. Once the 
triple (£ 1 , 5 , £ 3 ) has been determined by Rjj or Rjl, the idea is to ‘project’ the reflected 
wave along the associated integral curve; see Figure]^ (c), (d) for a picture of case z = 3. 


i) For z = 1, we compute Ru{U-,Ujf) = (£i,(5,£ 3 ) and let £1 commute with dg by 
The outgoing composite wave connects Up to t/+, where Ii{Sq){Up) = 


Lemma 


3.5 


d>i(£i)(17_) and U+ = l3{Sl + £3) o $ 2 (d) o $i(£i)(t/_). Hence, £0 = Rn{Up, U+) = 
(dg, S, (5q + £ 3 ). 

ii) For z = 3, we compute Ril{U-,Uj^) = (£i,(5,£ 3 ) and let £3 commute with dg by 
The outgoing composite wave connects U- to Uq, where Uq = I 3 {Sq) o 


Lemma 


3.5 


MS) O U{Sl + £i)([/_). Hence, £0 = Rii{U., Uq) = (S^ +8i,S, S^). 
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To prove (3.9), notice that for z = 1 we use Rix or Ru (i.e. di = h), while for z = 3 we use 
Rll or Rjl (i.e 63 = h). Hence, (3.9)2 is equivalent to 


a-9i{ei) + a+ 03 {e 3 ) 


a+h{ 6 i) 

a-h{ 63 ) 


if z = 1, 
if z = 3 . 


By (3.1) and (3.2) we have that for z = 1,3 


1 

2 


log 



a+h{di) 

a-h{53) 


Now, by Proposition 3.1 it suffices to notice that 


if z = 1, 
if z = 3 . 


£3 - ei = ^ log , 0-01 (ei) + a+ 03 (e 3 ) = ^ ^ 



Figure 6: Interaction of a 3-wave <53 with a composite wave Sq. (a), 
(c): the actual Riemann solvers, the Accurate case (a) and the Sim¬ 
plified one (c); (6), (d): the auxiliary pictures, the Accurate case (6) 
and the Simplified one (d). 


Finally, we verify the relations on the signs of the outgoing waves ei, £ 3 . We prove only case 
(5 > 0, since the other one is symmetric by replacing z = 1 with z = 3. Notice that for Rjx the 
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results collected in (3.10) have already been proved in [3] and we obtain the same interaction 
patterns of [U (5.5)]. 


When the Simplified solver is used and z = 3, by (3.9) £ 1,63 solve 


£3 — El — <53 ) 

a_£i + a+h{ez) = a-h{Ss). 


(3.11) 


Substituting the expression for £1 coming from the first equation of (3.11) into the second one, we 


obtain e 3 +kh{£ 3 ) = S 3 +h{ 53 ), where /c = a+/o_ > 1. Hence, we have that sgn £3 = sgn( 53 . Now, 
take (^3 < 0 and assume to use Rll to solve the Riemann problem at the point of interaction. 
The corresponding outgoing waves £^, £3 solve 


'3 “ — <^3 , 

a-h{el) + a+h{e'^) = a-h{ 63 ). 


(3.12) 


Since £^ > 0, then system (3.12) reduces to (3.11) and, by uniqueness, its solution coincides 
precisely with £i,£ 3 . Hence, (3.10) is valid. If <^3 > 0, instead, we have that h{ 53 ) = 83 and 
/i(£3) = £3, i.e. in this case it holds Rjl = Rn - This amounts to solve a linear system in £i,£3 
and we find £1 = — 835 / 2 . Hence, sgn£i = —sgn <53 = —sgmJ • sgn 53 , as wished. 

When i = 1, by (|3.9|) £i ,£3 solve 


£3 — £1 = —(5i, 

a-h{ei) + a+£3 = aj^h{8i). 


(3.13) 


Again, it is easy to prove that sgn£i = sgn5i. If > 0, then Rjj = Ru and system (3.13) 


is linear. Thus, we get £3 = 5i5/2 and sgnes = sgnJi = sgnJ • sgn5i, as wished. If, instead, 


(5i < 0, then the second formula in (3.13) becomes 


sinh £1 + k £3 = k sinh (5i, 


(3.14) 


where A: = a+/o_ > 1. By substituting the expression for £i obtained from the first equation of 
( [3T^ in ( [3T4l ), we get A :(£3 + (5i)+ sinh (£3 + 5i) = k (sinhJi + 5i). If we call r(x) := A:x + sinhx, 
then r(e 3 + 81 ) = k (sinh(5i + <5i) and 


r (£3 + (^ 1 ) - r((5i) = (A; - 1) sinhJi 


(3.15) 


Since T is a strictly increasing function and < 0, it follows £3 < 0, that is sgn £3 = sgn(ii = 
sgn J • sgn(5i. Therefore, the proposition is completely proved. □ 

Remark 3.8. Differently from in the Simplified Riemann solver the emerging error is not 
only on the u-component of the 0 -wave and the transmitted i-wave £i does not maintain the 
same strength 8 i of the incoming one. The latter is a key feature of the solver, that guarantees 


the decrease of the functional defined in (5.1) across any interaction. Indeed, here we take into 


aecount the possible appearance of a reflected j-wave £j (j i) that we attach to (5o in place of a 
standard non-physical wave as in this is possible because the states connected by £j and by 


the 5r. component lie on the same j-integral curve. See Remark 5.2 for more details. 
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4 Approximate solutions 


We use Proposition |3.7| to build up the piecewise-constant approximate solutions to (1.1) that 


are needed for the wave-front tracking scheme Pi- We first approximate the initial data ( |1.3D : 
for any z/ G N we take a sequence {Vg,Ug) of piecewise constant functions with a finite number 
of jumps such that, denoting 


1. TV (log(p^)) < TV (log(po)), TV «) < TV (u^); 


2. Iim2;^_oo«,0(x) = lim3,^_oo(uo,Uo)(x); 

3 - ll«,0 - (w,Wo)||li < l/i^- 


We introduce two strictly positive parameters: a = ay, that controls the size of rarefactions, 
and a threshold p = py, that determines which of the two Riemann solver is to be used and 
depends on the initial data. Here follows a description of the scheme that improves the algorithm 
of [Ij and adapts it to the current situation. 


(i) At time t = 0 we solve the Riemann problems at each point of jump of (Uq , Ao)(-, 0+) 

as follows: shocks are not modified while rarefactions are approximated by fans of waves, 
each of them having size less than a. More precisely, a rarefaction of size £ is approximated 
by iV = [e/cr] -|- 1 waves whose size is e/iV < a; we set their speeds to be equal to the 
characteristic speed of the state at the right. Then (u^, Ao)(-, t) is defined until some 

wave fronts interact; by slightly changing the speed of some waves we can assume that 
only two fronts interact at a time. 

(ii) When two wave fronts of the families 1 or 3 interact, we solve the Riemann problem at the 
interaction point. If one of the incoming waves is a rarefaction, after the interaction it is 
prolonged (if it still exists) as a single discontinuity with speed equal to the characteristic 
speed of the state at the right. If a new rarefaction is generated, we employ the Riemann 
solver described in step (i) and split the rarefaction into a fan of waves having size less 
than a. 


(in) When a wave front of family 1 or 3 with strength 5 interacts with one of the composite 
waves at a time t > 0, we proceed as follows: 


if |<^| > P, we use the Accurate solver introduced in Proposition 3.7 partitioning the 
possibly new rarefaction according to (i)] 


if Id I < p, we use the Simplified solver of Proposition 3.7 


5 Interactions 

In this section we analyze interactions between waves. We separately study interactions that 
involve one of the two composite waves and interactions between 3- and 1-waves occurring in 
one of the regions £, At,77. In particular, we focus on the interaction estimates for the former 
ones and we introduce a new functional F, different from that of jl], to estimate the possible 
increase of the total variation. 


13 






Consider i > 0 at which no interactions occur and (^ > 1 to be determined. Using indices 
m, r to refer to waves in the region C,M,TZ, respectively, we define L = + L"* + U, where 


i=l,3,5^>0 1=1,3,5^<0 

6^GC,M,n S^^C,M,TZ 


For > 0 (see Figure 


7) we introduce Q 


Q^ + Q'^ + Q\ where 


= (iF>| + iF^lCl) 1^31 + ’ 

53>0 <53<0 

5]^>0 53>0 

(5ieA1 6 ^^M. 

Q' = (a-;i„i + A-JICl) E l^il + E ■ 

5i>0 5i<0 

5 ]^ GlZ Si^TZ 


Moreover, we define and 

F = F^ + F’" + F^ + L°, (5.1) 

where = |r/o| + |Co|- Clearly, F can be seen as defined also by F = L + Q + F°. We also write 

l = L^ + L^ + lr =Y^ |<5*| + ^ ^ I<5,| = 2 TV (logp(t, •)) - hoi - iCol. 

i=l,3 i=l,3 i=l,3 

S^GM S^GTZ 


c 


M 


n 




Vo 




Figure 7: The parameters related to the approached 2-wave 

and to the regions of provenience of the approaching waves. 


Remark 5.1. The summation in (Q'^) is performed over the set of waves approaching the 

composite waves from the left (right, respectively) hut it does not include 3-shocks approaching 
Co (l-shocks approaching po, respectively). Moreover, the sum in Q'” includes neither 3-shocks 
approaching Co nor l-shocks approaching po. Indeed, the contributions given by these waves 
can be dropped from the interaction potential since the linear functional L decreases when they 
interact with a 0-wave. The functional F obtained in this way has proven to provide the best 
possible conditions on the parameters involved and, consequently, the largest ones on the initial 
data. Clearly, the choice of F is reflected in (2.3), where the total variation of the data can be 
taken larger in A4 than in the outer regions. 
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Remark 5.2. We will prove that the functional F decreases when we use the Riemann solvers 


introduced in Proposition 3.1. This property does not hold true with the solvers of As 


already mentioned in Remark 3.8. the key point is that in the Simplified solver the strength of 


the transmitted wave is not the same of the incoming one, while they coincide for the Pseudo 
Simplified solver ofm Indeed, consider an asymmetrical functional Fi adapted to the situation 
ofJHi, i.e. Fi does not include 3-shocks in the interaction potential. In the case of an interaction 
of a 3-shock with the 2-wave 62 > 0 solved by the Pseudo Simplified solver, we would get AL = 
-f^np| 72 ,o| > 0 and AQ = 0. Thus, Fi would increase. 

In the following we often assume that, for some fixed ttIq > 0, any interacting i-wave, i = 1, 3, 
with strength 5i satisfies 

\5i\<mo. (5.2) 

In particular, this bound is to be imposed only to shock waves, since we can control the strength 


of the rarefaction waves by ( 6 . 1 ) below. 


5.1 Interactions with the composite waves 

Here we collect all the estimates concerning the composite waves. 

Lemma 5.3 (Interaction estimates). Let i = 1,3. Consider the interaction of an i-wave 5i with 
a composite wave 60 = ((5g,(5, 5 q). Denote by Si the strength of the transmitted wave and by ej, 
j = 1,3, j 7 ^ i, the strength of the reflected one (even in the Simplified case, where it is attached 
to 5o). Then, when |(5j| > p it holds 


\ei - 5i\ = \ej\ <-\6i6\ and |eo-5o|=0: 


(5.3) 


while, when \5i\ < p it holds 


\^i hi\ — I £0 I — I £ J I ^ ^ 


' (J 

if 6 i < 0 and either (i = 1, 5 > 0) or (i = 3, 6 < 0), 


-|(5i5| otherwise, 


(5.4) 


where Co = Co{p) = sinh(/ 9)//9 > 1 is such that Co{p) —>• 1"'' for p ^ 0 '^. 


Proof. When |(5j| > p, i.e. when the Accurate solver is used, ( 5 . 3)2 is immediate and (|5.3[)i can 


be derived from (3.9) (case 01 = ^3 = h) following the same steps as in [3l Theorem 2]. 

When |(5j| < p, i.e. when the Simplihed solver is used, we analyze only the case 5 > 0 and 
refer to Figure]^ (c), (d). The equality |ej — (5j| = \ej\ in (5.4) is a consequence of (3.9)i, while 
ko ~ = kjl reflects our choice to attach the reflected wave to the composite one. To prove 


the inequality in (5.4), we distinguish cases according to the characteristic family and the sign 


of the interacting wave. If 5i > 0, then Ru = Ru for i = 1 and Rn = Ru for i = 3; moreover, 
it holds kjl = \5i5\/2 by (3.9). When i = 3 and the interacting wave has negative sign, as 


in Proposition 3.7 we have Rjl = Rll and the interaction estimate (5.4) follows exactly as in 


the Accurate case. Instead, when i = 1 and the interacting wave has negative sign, we have 
to pay more attention. Recall from Proposition |3.7| that we have £1 = £3 + and £3 < 0, 
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(5i < 0; moreover, ( |3.15 ) holds. By the Mean Value Theorem there exists some s such that 
r(e 3 + 5i) — r(hi) = r'(s)e 3 . Hence, we have 


and we deduce 


{k + l)|e 3 | < r'(s)|e 3 | = (A: - 1) sinh |hi| 


ksl < sinh|(5i| = ^smh|5i| < • 


□ 


We will discuss in Appendix [A| a refinement of estimate (5.4). Notice also that in the previous 
lemma the biggest effort is required to handle the estimates for shocks interacting with (5o and 
going towards the phase where a is smaller. In our case, these are precisely the shocks that hit 
Ai from the outside, i.e. 1-shocks interacting with ^md 3-shocks with ijq. 

Now, we are ready to give a first list of conditions to impose on the parameters and 

p in order that the functional F decreases at any interaction time. 

Proposition 5.4. Assume that at a time t > 0 a wave 5i, i = 1,3, interacts with one of the 
composite waves po or (q. Then, A.F{t) < 0 provided that 


^> 1 , 


aT(i + y 


2 S A'c < 1^1 . 


^^1 
2 - \p\ ' 


< + {Kl - 1)|C| , K^[l + ^)|C| < Kl\C\ + (A:^ - 1)|77| , (5.6) 


Co{p) < ^^min{A:^,A:^}. 


(5.5) 


(5.7) 


C 



n 


£ 


Vo 



(b) 


n 


Figure 8: Interactions of 1- and 3-waves with fo solved by means of 
the Accurate solver. Here the fronts carrying the composite waves 
are represented as a single line. 

Proof. Since the two cases give symmetric conditions, we only analyze interactions involving ^o; 
see Figure]^ We have 


£3 — ei = —di, 
£3 — Si = bs, 


|ei| - |(5i| = |£3| , if i = 1, 

kal - If^sl = -|£i| , if z = 3. 
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z = 1. If the interacting wave is a rarefaction, then by (5.3),(5.4) we have AL + A.L^ = l^al + 


ki| - l^il = 2|e3| < |5iC| and AQ = KJ^\eir]\ - K^|(5iry| - K^\ 6 iC\. Therefore, 


AF < Idil 


K 


(i + !|!) -A';)l,l + (i-A'j)ici 


which is nonpositive by (5.6)i. Instead, if the interacting wave is a shock, then by (5.3),(5.4) 


AL + AL“ = 


eki|+e|e3|-el<5i| = 2^1631 <^|-5iC| 

?kil + ksl ~ = (1 + OksI <^0(1 + ^) 


kiCI 


if ki| > p, 
if |(5i| < p, 


and AQ = —ia both cases. Consequently, 


AF < 


l-Kl 




kiCI ifki|>P, 

kiCI ifki|</=5, 


is nonpositive by (5.5) 1^2 and (5.7). 


i = 3. If the interacting wave is a rarefaction, then by the interaction estimates 

ksCI 


AL + AL^ = 


Ckil + kal - ksl = (^ - i)kil < (? -1)- 


kil + ksl - k 3 | = 0 

and AQ = -K^\ 63 C\. Then, 

■^-1 


if kal >P, 
if kal < P, 


AF 


< 


- 

2 ^ 


= -K^mi 


ksCI ifk3|>P, 

if ksl < P, 


is nonpositive if (5.5)3 holds. On the other hand, if the interacting wave is a shock, then 


AL + ALO = |ei| + ^|e 3 | - Cksl = -(? - l)ki| and 

i7™|ei7/| if|<53|>/3, 


AQ = < 


if ksl < P- 


Therefore, 


AF = 


-(e-l) + if™|r/|J |ei| if|53|>P, 

-(^-l)ki| if|53|<p. 


is nonpositive by (5.5)4. 


□ 
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Figure 9: Interactions of 1- and 3-waves. 


5.2 Interactions between waves of the same family 

In this subsection we analyze the interactions between 1- and 3-waves, see Figure 

For interactions of two i-waves, i = 1, 3, under the notation of Figure we shall make use of 
the identities [3]: 

{ —oi — /3i if i = 1 , 

h{ei) + hies) = h{ai) + h{l3i), i = l,3. (5.8) 

a3 + Ps if i = 3 , 

Lemma 5.5. For the interaction patterns in Figure^ the following holds. 


(i) Two interacting waves of different families cross each other without changing strengths, 
(a) Let ai, f3i be two interacting waves of the same family and ei, the outgoing waves. 

(a. a) If both incoming waves are shocks, then the outgoing wave of the same family is a 
shock and satisfies |ei| > max{|aj|, |/3i|}; the reflected wave is a rarefaction. 

(ii.b) If the incoming waves have different signs, then the reflected wave is a shock; both the 
amounts of shocks and rarefactions of the i-th family decrease across the interaction. 
Moreover for j i and ai < 0 < fli one has 

r*osn T _ 1 

lef < ciai) • mm{\a,\,m , c{z) = . (5.9) 

cosh z + 1 


The proof of Lemma 5.5 can be found in [T], where the function c is used in place of the 
damping coefficient d of [1]. Remark also that, by definition of the functionals, we need to 
distinguish between interactions taking place in Ai and in C or IZ. 


Proposition 5.6. Consider the interaction at time t > 0 of two waves of the same family 1 or 
3 and assume (5.2). Then, A.F{t) < 0 provided that 


1 <^< 


c(mo) ’ 


< 


KL 


+ Kl\C\<C-l, 


^-1 

ICI 

Kiln 


iL™ < ^^ 

^ - \ri\ ' 
TRiJlCI <e-l. 


(5.10) 

(5.11) 


Proof. First, we consider the interactions taking place in Ai, see Figure pR) Here, we only cover 
the case of interactions between two 3-waves 03 and fis giving rise to ei and £3 (the 1-waves 
case is analogous). 
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Figure 10: Interactions of 3-waves in Al. 


When both and are shocks, by Lemma 5.5 we have that Si is a rarefaction and we notice 
as in P3 Proposition 5.8] that 

AL+|ei|(^-l) = 0, (5.12) 

for any C ^ 1- Moreover, we have 

AQ = iL™|ei77| , 

-(C-i) + iL™H 


AF = 




and F is non-increasing by (5.10) 1 ^ 3 . On the other hand, when the two interacting waves are of 
different type, for example 0:3 < 0 < /Ss, as in PQ Proposition 5.8] one can prove that it holds 


AL + ^(e-l)|ei| <0 


(5.13) 


by condition (5.10)i. If S 3 is a rarefaction, then AQ = (leal — l/^al) |C| and F decreases 
by Lemma 5.5 if ea is a shock, then AQ = and, again, F decreases. Remark that 

the analysis of the interactions between 1 -waves requires symmetrically the condition iP™ < 

(C-1)/IC|. 

Next, we analyze the case of interactions taking place in TZ (similarly one proceeds to analyze 
those occurring in C). As for interactions between two 1-waves, it is easy to verify that F 
decreases with no need of other conditions than (5.10)i. 


£ 


M 


Vo 


^0 



Figure 11: Interactions of 3-waves in TZ. 


As for interactions between 3-waves (Figure [IT|), instead, we need to require (5.10)i^3 in order 
to have AF{t) < 0. Indeed, we have (|5.12[) when the interacting waves a 3 , P 3 are both shocks. 


while in the other two cases of interaction it still holds (5.13) under condition (5.11 )i. Also, if 
0 : 3 , /^a < 0 we have 


AQ = A;|eir?|+R^neiCl 


AF = 


- 1) + K^^\rj\ + Kl\C\ 


Nil > 
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while if, for example, as < 0 < /Sa we have 


AQ = Kl^\eiC\, 

AF<C\-{C-1) + KI\C\ 


l^il 


Consequently, F is non-increasing by (5.11 )i. 

Symmetrically, in C we get the condition K^\r]\ + K^(\C\ < ? — 1- 

5.3 Decreasing of the functional F and control of the variations 


□ 


In order that F decreases across any interaction, the various parameters in (5.2), (5.5), (5.6), 
(5.7), (5.10) and (5.11) are chosen in the following order. Given nio to be fixed later on, we 
choose in turn and finally p. Remark that in the following calculations we keep 

(almost) everywhere strict inequalities, since they are needed in the analysis on the control of 
the size of the composite waves (see Section 6.1). 


We notice that, for the choice of by (5. 5 ) 3,4 and (5.10) it must hold 


C -1 


< mm 


e-1 e-1 

\v\ ’ ICI 


(5.14) 


which is always satisfied since |f?|, |C| < 2. Moreover, by putting together the conditions obtained 

^ < KJ^\v\ + {Kl - l)|Cl < (e - 1) - Kl\C\ <(?-!)- Id • (5.15) 


in ( 5 . 5)3 with ( |5.6| )i we get necessarily 

n\H 


(e-i)(i + 

Hence, it follows 


K-i)(i + l|)l|l<({-i)-|C|, 


which is equivalent to 


1 + 


Id 


l-(l + |d/2)|d/2 




(5.16) 


provided that 1 — (1 -|- |d/2)|?/|/2 > 0. Analogously, from ( 5 . 5)4 and (5. 6)1 we get 

H 


1 + 




l-(l + |d/2)|d/2 

provided that 1 — (1 -|- |r 7 |/ 2 )|d /2 > 0. Therefore, it must hold 

Id Id 


1 -|- max 

under the stability condition 


l-(l + |d/2)|d/2’l-(l + ld/2)ld/2 




(5.17) 


|d\M 1 


min<!l-(l + ^j^,l-(l + 


2/2 


> 0 , 


which is equivalent to (2.1). Then, by (2.2), (5.5)i and (5.17) we obtain the condition 


i+'R(id,id)<e< 


1 


c{mo) 


Summarizing, the choice of the parameters proceeds as follows. Let d|, iCl satisfy ( 2 . 1 ). 
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Recalling (5.10)i, we fix mo such that 

c(mo) < 


1 


i+n\v\,\c\) 


(5.18) 


We will prove in Section [6.2| that this is possible under suitable conditions on the initial 
data. Notice that c is a strictly increasing function of mo and then it is invertible. 


Then, we choose ^ in the non-empty interval given by 




1 


c{mo) 


(5.19) 


We choose iT™, such that 

< RT™ < min 


^ ^ ^ < min 


e -1 (e-i)-ici 1 _ (c-i)-ici 

\v\ ’(r+|C|/2)|r/|/ (l + |C|/2)|r?|’ 

e-1 (e-i)-hi 1 ^ (g-i)-i^i 
Id ’(i + id/2)id/ (i + id/2)Kr 


(5.20) 


(5.21) 


This is possible since these two intervals are non-empty by ( |5.19 ) and ( |5.16 ). Thus, (5.5> 3^4 
and ( 5 . 10 ) 2,3 follow and it holds 

K,”(i + lfl)i„i < K -1) - ici, A';”(i + h!)ici <({-!)- Ill. 


(5.22) 


By (5.22), we choose K^,K^ that satisfy 

a”(i + 1|)i»i<a;h <«-!)-Id, 
N”(i + y)|‘;Isa'ICI<({-i)-I"?I; 

then, we can take and such that 


1 < < 1 + 


(e-i)-ICI-i^,^ 
“ICI ■■■ 


1 .... .1 , (e-D-M-^c^lCI 
1 < < 1 +--, 


(5.23) 

(5.24) 

(5.25) 

(5.26) 


from which (5.6) and (5.11) follow. 


• Finally, we choose p that satisfies (5.7). 

In the following proposition we collect the results obtained so far. 


Proposition 5.7 (Local decreasing). Let mo > 0 satisfy (5.18) and consider the interaction of 


two waves at time t that satisfy ( |5.2| ). Moreover, assume that d Krjf ’ P satisfy (5.19)- 
(5.26) and (5.7). Then, 

AF{t)<0. (5.27) 
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Now, we can prove the global in time decreasing of the functional F. 


Proposition 5.8 (Global decreasing). Let rrio > 0 satisfy (5.18) and choose parameters 
^ Proposition 


5.1. Moreover, assume that 


L^(0) + c{mo)L^{0) + L^{0) < moc{mo 


(5.28) 


and that the approximate solution is defined in [0, T]. Then, F{0) < mo, ^F{t) < 0 for every 


t G (0,T] and (5.2) is satisfied. 

Proof. For convenience, we use notation L'ffi and Lff to indicate the partial sums in L™' due to 


i-rarefaction waves [iR) and i-shock waves {iS), respectively. By ( 5 . 5 ) 3,4 we have 

^*"(0) = L’”(0) + g™(0) < 

< L’£(0) + LfslO) (1 + K”K|) + iS,(0) + LTn(O) (l + 

< L’,”s(0) + (LTMO) + LS(0) + (t-TMo) < {£”(0) ■ 


< 


Moreover, from (5.11) it follows 


F^'^(O) < L^’’'(0) (1 + K^/\ri\ + iFjnCi) < . 


Then, 


F(0) = F^(0) + F”^(0) + F’’(0) < fL\0) + eL™(0) + . (5.29) 

Now, for a fixed t <T, suppose by induction that F{t) < mo and AF(r) < 0 for every 0 < r < t, 


interaction time. Then, by Proposition 5.7 we have AF’(t) < 0. This implies 

Fit) < F{0) < fl^O) + |L™(0) + fL'^iO). 


By (5.29) and (5.28) the size |(5j| of a shock (i = 1,3) at time t satisfies 


|<5*| < yT(t) < CL^iO) + L™(0) + eT'-(O) < 


1 


c(mo) 


Z^(0) + Z”^(0) + 


c(mo) 


L^(0) < mo. 


Hence, (5.2) is satisfied and the proof is complete. 


□ 


6 The convergence and consistency of the algorithm 


In this section we finally conclude the proof of Theorem 2.1 on the convergence and consistency 
of the front tracking algorithm. 

In order to be well-defined, the algorithm must satisfy three main requirements: i) the size of 
rarefaction waves must remain small; ii) the total number of wave fronts and interactions must 
be finite; hi) the total size of the composite waves must vanish as the approximation parameter 
n tends to -|-oo. The first one is accomplished as in [H Lemma 6.1] and, in particular, the size 
e of any rarefaction wave is bounded by 


0 < e < cj ^1 -I- ^ max {|r/|, ICI}^ <2a. 


( 6 . 1 ) 


The second requirement can be proved as in [H Lemma 6.2 and Proposition 6.3]; while the next 
section is devoted to the proof of iii), i.e. of the consistency of the algorithm. 
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6.1 Control of the total size of the composite waves 

Our wave front-tracking scheme exploits the concept of generation order of a wave to prove that 
the size of the error attached to the two phase waves tends to zero. We introduce such generation 
order as in [3] for 1 - and 3-waves, while for the composite waves we proceed in the following 
way. Consider an interaction with 6 q (either rjo or (^o)- we assign order 1 to the component <5 
(which never changes) and order A:.y -|- 1 to the outgoing j-th component when the interacting 
wave is an z-wave 7 of size < p, i, j = 1,3, j ^ i, while we keep the order of the other component 
unchanged. In other words, denote by (fc^i, 1,the triple made of the orders of the three 
components and consider the interaction with a wave 7 solved by the Simplified solver. Then, 
for the outgoing 0 -wave sq = (ej, S, ej]) we have 


I (/c.y-hl,l,/c 53 ) if 7 is of family 3 , 

I (k^i ,l,k^ + 1 ) if 7 is of family 1 . 

Mimicking what has already been done in mm, for every k = 1 , 2 ,... we define the functionals 
Lk, Qk and Fk simply by referring L,Q,F to waves with order k. Moreover, we define Lk = 
Ylh>k^h ^'^d Fk = Ylh>k^h- In detail, for A: > 1 and for as in Proposition 5.7 we 

define Lk = Ll + + LI, where 


L 


£,m,r _ 


E 1*1+{ E 




i=l,3, Sj_>0, kg,=k 
6 ^£C,M,n 


z=l,3, 5^<0, kg, =k 
Sj_GC,M,n 


and Qk = Qk + Qk + Qk, where 


gi= (K^| 7 |+i^J|Cl) E 1 - 531 +E l^ 3 r/|, 


<53 ez: 


^3<0- '= 53 ='“ 
■53 e£ 


Q^ = K^ ^ \6,p\+K^ l^sCI, 


« 1 > 0 . kg,^=k 


^2 C.A/I 


Ql = iK^^\v\ + Kl 


Y \Si\+CKl Y 

5x>0,A;^^—fc fc 

<5^ G'R- G'R- 


Moreover, we define 


where 


Fk — Lk + Qk + ) 

7 = E (i^wi + i^ioi) ' 

Tk<t 


( 6 . 2 ) 

(6.3) 


with Tk in (6.3) denoting the interaction times when a small reflected wave of order k is born 


and attached to one of the composite waves. As a consequence, only the times where the 


Simplified solver is used give positive summands in (6.3), since when the Accurate solver is used 
we have \Apo\ -|- |ACo| = 0 . 

By Ik we denote the set of times in which two waves of the same family of order at most k 
interact with each other, while by Jk we denote the set of times in which a wave of order k hits 
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one of the two composite waves. Moreover, let Tk = /fc U Here, we will prove the analogous 
of [H Proposition 6.1]. In particular, we define: 


/i = max 


2^-1' 2Kl - r 1 + 2K^ ’ 1 + 2K^ 


l + 1+K^\C\ 


l + (it'^|r?|+/i'^^|C|) l + {K-^\v\+Kl\C\) Co 


Co 


(6.4) 


? ’ e ’ ’ e(2K^"-l) ]■ 

We have that fi < 1 hy the conditions required in Proposition |5.7[ 

Proposition 6.1. Let mo,^, satisfy the assumptions of Proposition 


5.1 and assume that 


(5.2) holds, that is, |<5i| < mo for the size of every wave. Then, the following holds for r G Th, 
h > 1: 


^Fh < 0 ) ^Fhj^i > 0 , 
APfc = 0 if k > h + 2 . 


Moreover, 


h-l 


[AFh+i]+ < MAFh]- -Y^AFi 


£=1 


(6.5) 

( 6 . 6 ) 


(6.7) 


Remark 6.2. Proposition 6.1 let us improve Proposition\5. 1\ Indeed, reealling that Fh = Ih^Jh, 


Proposition 6C_ implies, for t £ F, 

h-l 


AF = Y^Fe- + [AF/^+l]+ < -(1 - L)[^Fh]- <0, 


£=1 


while for r G Jh, being [APr]+ = 0, it gives 

AF = -[AF;,]_ + [AF;,+i]+ < -{I - yi)[AFh\- < 0. 


Then, estimate (6.7) specifies the deerease of the funetional F and improves (5.27). 


Proof of Proposition\6.1\ If A: > /i + 2, no wave of order k is involved and then (6.6) holds. To 


prove (6.5) and (6.7), we distinguish between two cases. 


r G I/i (Interactions between waves of 1-, 3-family). 

Clearly the F^s do not vary when a 1-wave interacts with a 3-wave. Then we consider 
interactions of waves of the same family occurring in one of the three distinct regions C,Ai,IZ. 


Since t G Ih, then ALh+i = AT 


£,m,r 

h +1 


> 0 and 


AQ)]Yi — AT)(Yi for interactions in M, 


0 < AQh+i = 


( 6 . 8 ) 


AQ'4i < (K. 


r\ 


+ K, 


A r^,r 
^^h +1 


for interactions in C, TZ. 


Also, AFh = ALh + AQh < 0, since both terms in the sum are negative or zero. This proves 

d^. 
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By (5.12) and (5.13) (see also [H (6.10)]), we have that 


1 / ^ ^ 

[AL,+i]+ < -([ALh]--Y,AL 


e=i 


(6.9) 


From (|6.9|), (| 6 . 8 |) and (|6.4|) we deduce that 

{ 


0 < AFh+i < < 


1 + max <J K^\r]\,K^ 


l + K^/\ 


+ K] 


h-1 


i,r 


|C||^ [AL/i+i]+ for interactions in M , 
[ALh+i]+ for interactions in £, 7^, 


( 6 . 10 ) 


<^l([ALh]--Y,^L^ 

We now claim that 


i=i 


h-1 


[^Qh\- — ^ > 0 • 


( 6 . 11 ) 


1=1 


To prove (6.11), we only have to analyze the cases when AQi^ > 0 for an £ < /i — 1. This can 
occur for interactions between waves a* < 0 < /3i with ctj of order i and /3i of order h, giving 
rise to waves Ej of different sign, i,j = 1, 3, z 7 ^ j. More precisely, by Lemma 5.5 we have 


h -1 


[AQ;,]_-^AQ, 


= 


1=1 


-ksl + m) {k^,\v\ + ^c'Kl) + eA'^l«3r/| 

iL-(-|ei| + |/3i|) |r?| 

(-1^3!+ 1/331) Id 

(-kil + |/3i|) (iLdd +^c"l^l) +^Kl\aiC\ 


for interactions in £ , 
for interactions in Ai , 
for interactions in Ai , 

for interactions in TZ , 


which is always a nonnegative quantity. This proves (6.11). Therefore, for r G Ih, estimate (6.7) 
follows from ( 6 . 10 ) and ( 6 . 11 ). 


T ^ Jfi (Interactions with the composite waves). 


Here we focus only on interactions involving since the other case gives symmetric condi¬ 


tions. Since no wave of order < h — 1 interact, then (6.7) reduces to 

[AF/i+i]+ < ^[AFh]- . 


( 6 . 12 ) 


To prove (6.12), assume that a 1-wave hi of order h interacts with (^oi see Figure 12 for the 
Accurate case. If hi is a rarefaction, then by Lemma 5.3 we have 

m ^ Idd 


AFh = ALh + AQh + ALl < 


<|d| 


Id 


+ iL-|eid-iFdhid-i^d'^iCI< 


Id 
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Figure 12: Interactions of a 1-wave di with solved by the Accurate 
solver. 


which is nonpositive by (5.6)i and, in particular, implies [AFX]_ > (2iF^ —l)|(5iC|/2. Notice that 
in the formula above the first inequality is due to the possible presence of a wave of generation 
order h in AL°. As a consequence, we get (6.12) by (6.4): 


[AF/i+i]+ = Lh+i = leal < < 


2 Kl - 1 


[AF,]_ 


If, instead, Ji is a shock, then 

/\Fh = i\e:,\-Kl^\6iC\+/\Ll<{ 




dCo-2K, 


M 

2 

|5iCI 


if li^il > P, 
if |(5i| < p. 


is nonpositive by (5.5)2 and [AF/i]_ > ^{2K^ — l)|(IiC|/2. Hence, by Lemma 5.3 we get 

1 


[AF/i+i]+ = < 


ih+i = ?k3l < 


ALO+1 = leal < ^l^iCI < 


a, 

2 


2a:^ - 


3;[AH/.]- 


if l^^il > P, 


Co 


-[AFh]- if|5i|<p. 


^{ 2 Kl - 1) 

On the other hand, let us consider the interaction with a wave 63 of order h belonging to 

we have AFX = ALh + AQh + AL^ < 


family 3. We first analyze the case 63 > 0; by Lemma 


5.3 


-[1-1- 2iL™]|ei| < 0. Thus, by (6.4) we get (6.12): 

Lh+i = ■Ckil < 

AT? , 1 = |ei| < 


[AF/i+i]+ = < 


1 + 2 K^ 


1 


1 + 2 K^ 


[AFh]- if|(Ia|>/0, 

[AT,,]_ if|(Ia|</0. 


In the other case, i.e. when 83 is a shock, we have AFh A ~'^|ei| < 0. Hence, 

i + iLr'-i 


[AF/i+i]+ = < 


Lh+l + Qh+l ~ (f l^ll 


< 




1 , 


AL^,+i = |ei| < -[AF,]_ 


'''^[AF;,]_ if|(Ia|>P, 
if |(5a| < p- 


Then, (6.12) is completely proved. Finally, we notice that in all above cases for t ^ Jh formula 
holds. □ 
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The remaining analysis aims at proving that 


j>k 


(6.13) 


for any k > 2 and for any t; it is carried out as in [H Propositions 6.3 and 6.4] and follows 
from Proposition 6.1 In particular, formula (6.13) is needed to prove that the total size of the 
composite waves vanishes as —)• oo. We conclude the section by determining parameters p and 
u as in [Ij. Fix a > 0 such that cr = (T;y—)-0asi ^—>-00 and estimate the total number of waves 
of order < k. Then, the total size of the composite waves is less (or equal) than 


Lk{t) + ^ (|Ar/o| + |ACo|) (t/i) 

hKk 

< • -^i(O) + C'o(/9)^ (|??| + Id) [number of fronts of order < k] 

< ■ ruo + Co{p)^ (|?7| + Id) [number of fronts of order < k] , 

which is < Ijv by choosing k sufficiently large to have the first term < l/{2u) and, then, 
p = py{mo) small enough to have the second term also < l/[2u). 


6.2 End of the Proof of Theorem 2.1 and a comparison 


In this last section we accomplish the proof of Theorem 2.1 and compare the result we obtain 
with that proved in mu- 


End of the Proof of Theorem 2.1. It only remains to reinterpret the choice of the parameter 


roo in terms of the assumption (2.3) on the initial data. Notice that we can approximate 


the initial datum (already satisfying 1., 2. and 3. of Section]^ in such a way that the jump 
[{p£,Ui, Xg), {pm,Um, Xm)) at the interface x = a is substituted by a jump consisting of the 2- 
wave separating {pi, U(, Xi) and {pe, u^, Xm) at x = a and by the solution to the newly appeared 
Riemann problem at x = with states {pi,U£, Xm) and {pm,Um, Xm)- Analogously, we can 
proceed for a jump [{pm,Um, Xm), {Pr,Ur, Xr)) at X = b. This is possible because p,u remain 
constant across a phase wave. Thus, we can relate hypothesis (2.3) to ( 5.28| ) by including in Z™ 
{U, respectively) the total variation of po and Uo at the interface and as in [U (3.12)] we can 
prove that 

T'(0)<^Tv(log(po),-) , 

2 x<a J 

(log(po), 


u 


‘(0) < - TV , _ 

2 a<x<() y 

1 


Up \ 
O-m J 


(6.14) 


T^(0)<^Tv(log(pd,- 

2 x>b y CLj, 


Now, by (5.19), (5.28) and (6.14) we have to look for an mp satisfying 


^(ld,ld)< 


1 


c{mo) 


- 1 = 


cosh rUo — 1 


=: w{mp) 


(6.15) 
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and 


TV 

x< 


I (log(po), + c{mo) TV (log{po), + TV (\og{po 

O \ 0,1 J a<x<6 y Qm J x:>6 y 


Ur 


), — < 2moc{mo) =: z{mo) ■ 

dr J 

(6.16) 

Notice that w{mo) is strictly decreasing from M_|_ to M+, while z(mo) is strictly increasing on 
the same sets. We can now define 

/C(r) = z , r G (0,+oo), (6-17) 

which can be written explicitly as 

= = + + ■ ( 618 ) 

Hence, if the assumption (2.3) holds, it is easy to prove that one can choose ruo such that (6.15), 
(6.16) hold. Finally, in order to pass to the limit and prove the convergence to a weak solution, 
one can proceed as in [9]. Theorem 2.1 is, therefore, completely proved. □ 

Remark 6.3. Notice that a slight improvement of Theorem \2.1\ follows from the use of the 
Riemann invariants 

z = u — a log n , w = u + a log v , 

where a = a^, a = am and a = in the regions C, M. and TZ, respectively. Indeed, recalling 
Definition 3.2. one easily finds that the solution to the Riemann problem with U- = (n_,M_, A) 
and U+ = (n+,u+,A) satisfies 

ki| + lesi < ^ (|'W)+ - w-\ + | 2 + - z-\) < ^1 log(p+) - log(p_)| + ^\u+ - u-\ , 

(with obvious notation) and the second inequality is possibly strict (for instance if the solution 
to the Riemann problem is a single rarefaction). 

Hence, the right sides of (6.14) could be replaced by TVa;<fl (wo, Zo), 4 ^ '~INa<x<6 {wo, ^o) 
4 ^ TVx>i (iCo, 2 (o); respeetively, and (6.16) could be given in terms of these quantities, 
leading to a weaker assumption on the initial data. 


Now, we make a comparison between Theorem 2.1 and the main result in [3j, which was 
proved to be equivalent to Theorem 3.1 of [5]. First, notice that condition (3.6) of [5l Theorem 
3.1] can be written as \p\ + jd < 1/2, when applied to the current problem. Then, it implies 
(2.1), since 




(i+¥' 

hi 

\ V 2 ^ 

' 2 ’ 

V 2 y 

' 2 / 


+ ICI 


<4’ 


i.e. the domain D contains entirely that of [5]. 


Next, we claim that H{\r]\, |C|) < Ivl + |CI when \p\ + jCj < 1/2. Indeed, by (2.1) we have that 

ICI 


is equivalent to 


l-(l + |C|/2)H/2 

ICI 


< hi + ICI 


+ ICI)(i +y) <2, 
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which holds true by (2.1). Similarly, the inequality 


l-(l + |r?|/2)|C|/2 


< \v\ + ICI, 


is equivalent to (|r/| + |C|)(1 + |^|/2) < 2, which holds true for the same reason. This proves the 
claim. 


On the other hand, condition (3.7) of [5l Theorem 3.1] here becomes by (1.4) 


TV (log(po), —Uo^ < H{\ri\ + Id) 
\ / 


(6.19) 


where the function H(r) is only defined for r < 1/2 by 


H{r) = 2(1 — 2r)k ^(r), k{mo) = 


1 - y/d{mo) 

2 - y/d{mo) 


( 6 . 20 ) 


Here above, d{mo) is the damping coefficient introduced in [U Lemma 5.6]. From [T] we already 
know that JC > H in the common range \r]\ + |d < 1/2. Thus, (2.3) improves ( |6.19 ), since the 

and 


left-hand side of (6.19) is bigger than the left-hand side of 

/C(77(|r/|,|Cl))>/C(|d + ld)>^(kl + ICI). 

Consequently, we obtain enhanced conditions on the initial data in comparison with mis] , even 
though the latter results apply to a wider class of Aq. 

Remark 6.4. As already mentioned, in absence of one of two phase-waves, i.e. for example 


r] = 0, hypothesis |21 (2.3)], as well as (2.4), can he improved by (2.5). Indeed, consider a 


modified functional F = F™ -|- F’’ -|- L^, where F™’^,F^ are the same of (5.1) with r/ = 0 and 
< b}. An 

c{mo) < 


M = {{x,t) : X < b}. Arguing as in the proof of Proposition 5.8. if we assume that 

1 


and c(mo)L”^(0)-b L^(0) < moc(mo) 
then F decreases. This last inequality holds true if we take 

c{mo) TV (log{po), —) -b TV (log{po), —) < z{mo) , 
X<CB Ojyyi J B Qj'f' J 


which is implied by (2.5). 


A On the estimate (5.4) 


In this short appendix we consider estimate (5.4), in the case of interaction of a 1-shock with a 


composite wave Jq. The question is whether we can snbstitute Co with 1. The answer is only 


partially positive as we show in the following lemma, which improves (5.4). For simplicity, we 
focns on the case (5 > 0. 
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Lemma A.l. Assume that a 1-shock (^i interacts with Sq, when 5 > 0 and |5i| < p. Then, 


ko - <^o| = ksl < */0 < (5 < \/5 - 1, 


where the bound on 6 in (A.l) is sharp. We also have 

\eo-So\ = \e3\<^^^6\6,\, 


(A.l) 


(A.2) 


for a suitable function 0 such that 0(5, z) < 1 for 6 < 2/3, Q{6, z) > 1 on a right neighborhood 
of z = 0 for S > 2/3 and lim^-^+oo 0(5, z) = 0. 


Remark that conditions (A.l) and (A.2) overlap when 5 < 2/3. 

Proof of Lemma \A.1\ First, we claim that for every fixed 5 > \/5 — 1 there exist r = r(5) and 
R = R(S) such that 


sup 

|5i|<r 


ko — <5o| 
5|5 i|/2 


> 1 , 


sup 

|<Si|>ii 


|eo — (5o| 
<5|5 i|/2 


< 1 . 


(A.3) 


Recalling that in this case ei = £3 + 5i and 63, 5i < 0, equation ( |3.14 ) can be written as 

A:|e3| + sinh(|e3| + |5i|) - /csinh |5i| = 0 , 

i.e. r(A:, x,g) = 0 where x = |5i|, y = |e3|, k = (2 + 5)/(2 — 5), r(k,x,y) = ky -h sinh(x + y) — 
/ssinhx; this equation implicitly defines y = y(x;k). Then, ( |A.l ) for a fixed 5 is equivalent to 
y < and, since Ty = k + cosh(x + y) > 0, this is true iff T{k, x, > 0 for all x > 0. R 

holds 


V [k,x, 


k-l 
k + V 


x\ = {k 


k-l 


+ 


2k 


k -\-l fc + 


- -fcj x + ^ 


h=l L 


2k 

fc + 1 


2h+l 


- k 


,. 2 / 1+1 


(2h + 1)! 


=: ^ ak,h 


,. 2 / 1+1 


h=l 


(2/1 + 1)! ’ 


as well as 


, A: — 1 \ 3 ,3. k{k — l){k‘^ — 4A: — 1) 3 3 

T [k,x, -——X ) = ttkpx^ + o{x'^) =-——- -x^ + o{x'^). 


k + 1 


6(/s + l)3 


We have < 0 iff /c > 2 + \/5; hence, f(A:, x, |^x) < 0 in a right neighborhood of x = 0. 
Since k > 2 + y/5 iff 5 > VE — 1, (A.3)i follows. Moreover, for /c > 1 fixed, there exists h{k) 
such that ak,h > 0 for /i > h{k)] this implies lim^j^+oo r(A:, x, ^f^x) = +00, hence (A.3)2. This 
proves the claim. As a consequence, for 5 > ^/E — 1 it holds 


sup 

I++0 


ko ~ <^ol 

<5|5i|/2 


> 1 , 


and the estimate on the left in (A.l) fails. 


0n the contrary, when /c < 2 + ^/E (i.e. 5 < VE — 1), there holds 0 < 0^,1 < ak,h for h > 1, 
k-l , 
fc+i' 


hence T{k,x, tttx) > 0 for all x > 0 and (A.l) follows 
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Now, we prove (A.2); by ( 3.14[ ) and the Mean Value Theorem, there exists s G]e 3 + 
such that r(e 3 + (5i) — r(hi) = r'(s)e 3 . Since s < hi < 0 we get r'(s) > r'(hi) = /c + cosh hi so 


that (A; + cosh I hi I) • |e 3 | < {k — l)sinh|hi|. Then, (A.2) holds for all h > 0, since 
A:-l 


ks < 


k + cosh |hi| 


sinh |hi| = 


A: + 1 sinh I hi I 
A: + cosh|hi| |hi| 


^|hhi|=: 0 (h,|hi|).^|hhi|. 


Finally, let us prove the properties of 0. For simplicity, call z = |hi|; then 0(h, z) < 1 iff 
0 < (A: + coshz) — (A: + 1 ) sinhz/z, which is equivalent to 


[(2h + 1) — (fc + 1)] 


^2h 


h=i 


(2h + 1)! 


> 0 . 


This last inequality is verihed for every 2>0iff2/i + l>A: + l for every h > 1, i.e. iff A: < 2. 
It is easy to check that A: = 2 is equivalent to h = 2/3. This implies that 0(h, z) < 1 for every 
2 ; > 0 if h < 2/3, while 0(6, z) > 1 on a right neighborhood of 0 if h > 2/3. By construction 
0(h, hi) < C'o(hi) < Co(p), with lim^^^oo ©(k <^i) = 0; while lim 5 ^_,.oo C'o(hi) = + 00 . Hence, 
(A.2) is far better than (5.4) especially when h > 2/3. □ 
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